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ABSTRACT 

We report the detection of absorption lines by the reactive ions OH + , H 2 + and H 3 + along the line of sight to the submillimeter continuum source 
G10.6-0.4 (W31C). We used the Herschel HIFI instrument in dual beam switch mode to observe the ground state rotational transitions of OH + 
at 971 GHz, H 2 + at 1115 and 607 GHz, and H 3 + at 984 GHz. The resultant spectra show deep absorption over a broad velocity range that 
originates in the interstellar matter along the line of sight to G10.6-0.4 as well as in the molecular gas directly associated with that source. The 
OH + spectrum reaches saturation over most velocities corresponding to the foreground gas, while the opacity of the H 2 + lines remains lower than 
1 in the same velocity range, and the H30 + line shows only weak absorption. For Local Standard of Rest (LSR) velocities between 7 and 50 kms~' 
we estimate total column densities of iV(OH + ) > 2.5 x 10 14 cm" 2 , N(H 2 + ) ~ 6 x 10 13 cirr 2 and N(H 3 + ) ~ 4.0 x 10 13 cm 4 . These detections 
confirm the role of + and OH + in initiating the oxygen chemistry in diffuse molecular gas and strengthen our understanding of the gas phase 
production of water. The high ratio of the OH + by the H 2 + column density implies that these species predominantly trace low-density gas with a 
small fraction of hydrogen in molecular form. 
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1. Introduction 

Oxygen chemistry in the interstellar medium is initiated by 
cosmic -ray or X-ray ionization of hydrogen. The resulting H + 
and HJ transfer their charge to form the + and OH + ions, both 
of which react rapidly with H2 molecules. Subsequent H-atom 
abstraction reactions proceed to the terminal ion H30 + with high 
efficiency, unless the electron fraction is unusually high. H30 + 
does not react with H2 or any other abundant neutral species, 
but is destroyed by dissociative recombination with electrons. 
Recombination produces O, OH, and H2O. This sequence is 
the main gas-phase chemical route to water in the interstellar 
medium. Until now, this classical scheme of the oxygen chem- 
istry (Herbst & Klemperer 1973) has been tested mainly through 
observations of OH and H2O, whose abundances are affected 
by many other processes as well, and by a few observations of 
H30 + . With the HIFI instrument on the Herschel satellite it is 
now possible to observe the strongest transitions of light hy- 
drides, including the reactive molecular ions OH + and H20 + that 
are thought to play crucial intermediate roles in the formation of 
OH and H 2 

Previous attempts to observe interstellar H20 + and OH + 
from optical spectroscopy are very limited. Smith et al. (1984) 
report an upper limit of 10 12 crrT 2 for H20 + based on optical 
absorption spectroscopy. No observation of the OH + transition 
at 358.4 nm has been reported until very recently (Krelowski 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant particiation from NASA. 



et al., private communication). Submillimeter telescopes have 
been more successful in detecting these ions. With the Atacama 
Pathfinder Experiment (APEX) telescope Wyrowski et al. 
(2010) detected two hyperfine-structure components of the OH + 
N - 1-07 = 0-1 transition near 909.2 GHz towards the strong 
continuum source Sgr B2(M) and concluded that OH + plays a 
major role in the chemistry of the diffuse molecular gas. Using 
HIFI on Herschel, Ossenkopf et al. (2010) report the detection 
of o-H20 + absorption in the three massive star-forming regions 
NGC 6334, Sgr B2, and DR 21. 



We present observations toward the massive star-forming 
region G 10.6-0.4 (W 31C) located at RA = 18M0m28.70,s, 
Dec = -19°55'50s (J2000). The submillimeter-wave contin- 
uum emission is associated with a cluster of ultra-compact H 11 
regions located in the so-called "30 km s _1 arm" at a kine- 
matic distance of 4.8 kpc (Fish et al. 2003, Corbel & Eikenberry 
2004). The molecular cloud surrounding the star-forming com- 
plex shows evidence of continuing accretion onto the embedded 
stars (Keto & Wood 2006; Sollins & Ho 2005). As discussed by 
Neufeld et al. (2010), line absorption by foreground gas has been 
detected from [OI] (Keene et al. 1999), I3 CH + (Falgarone et al. 
2005), HCO+, HCN, HNC, CN (Godard et al. 2010), and atomic 
hydrogen (Fish et al. 2003). While the gas associated directly 
with G10. 6-0.4 is detected in the velocity range -5 to +5 km 
s , the foreground gas along the line of sight is confined to the 



velocity interval from ~ 5 to 45 km s 
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Table 1. Transition spectroscopic parameters 



Transition 


Frequency 


Error 


E, 


A 


Ref 




MHz 


MHz 


cnT 1 


10-V 1 




OH + N =1-0 












2,5/2- 1,3/2 


971803.8 


1.5 


0.0 


1.82 


1 


2,3/2-1,1/2 


971805.3 


1.5 


0.0 


1.52 


1 


2,3/2-1,3/2 


971919.2 


1.0 


0.0 


0.30 


1 



o-H 2 0+l u -0 ,o 
3/2,3/2-1/2,1/2 

3/2,1/2- 1/2,1/2 
3/2,5/2- 1/2,3/2 
3/2,3/2-1/2,3/2 
3/2,1/2-1/2,3/2 



1115122.0 
1115158.0 
1115175.8 
1115235.6 
1115271.6 



10 
10 
10 
10 
10 



0.0 
0.0 
0.0 
0.0 
0.0 



1.71 
2.75 
3.10 
1.39 
0.35 



p-H 2 + l li0 -lo,i 

3/2,3/2-3/2,3/2 


607207.0 


20 


20.9 


0.60 


2 


H,0 + 

- lo + 


984711.9 


0.3 


5.1 


2.3 


3 



1, Muller et al. 2005 & CDMS; 2, Strahan et al. 1986, Ossenkopf et al. 
2010; 3 Yu et al. 2009 & JPL catalog 



2. Spectroscopy 

Table [TJ summarizes the spectroscopic parameters used for the 
searched transitions. The OH + ion has a 3 2~ electronic ground 
state. We used the data from Beekoy et al. (1985) with re- 
ported uncertainties for the line frequencies better than about 
1.5 MHz. Muller et al (2005) have published a global fit us- 
ing these data which is available in the Cologne Data base for 
Molecular Spectroscopy (CDMS). The ground state dipole mo- 
ment obtained from theoretical calculations is 2.3 D (Werner et 
al. 1983, Cheng et al. 2007). 

The spectroscopic properties of H20 + are extensively dis- 
cussed in the accompanying letter (Ossenkopf et al. 2010). Its 
dipole moment is very similar to that of OH + , as recently cal- 
culated by Wu et al. (2004), who found the ground state dipole 
moment to be 0.933, mistakenly printed in atomic units for a 
value of 2.370 Debye. This molecular ion has ortho and para 
symmetry states, which have the opposite parity to the water va- 
por molecule, i.e. the lowest state is the ortho state. Therefore 
the ortho/para ratio should be at least 3:1. 

H30 + has been the subject of many spectroscopic investi- 
gations, because for the frequency of the inversion transition of 
this molecular ion is comparable to the rotational constant. The 
most recent analysis has been presented by Yu et al. (2009) and 
is available in the JPL spectral line catalog. We observed the 
lowest transition from the E symmetry state, which is ~ 5 ctrT 1 
above the ground state. 

3. Observations 

The observations were carried out by the HIFI instrument on 
board Herschel (Pilbratt et al. 2010; de Graauw et al. 2010), 
in the double beam switching (DBS) mode. We used the Wide 
Band Spectrometer (WBS) with 4GHz bandwidth and 1.1 MHz 
effective spectral resolution. To help separate the signals from 
the lower and upper sidebands, three integrations were carried 
out at each frequency, with the local oscillator (LO) shifted by 
lSkms" 1 . 

The data were first processed with HIPE dOtt et alJ2010h and 
subsequently exported to CLASS. We checked that the target 
lines did not suffer from any contamination by lines from the 
other sideband and that the absorption features were not con- 
taminated by emission in the reference beam. The two polar- 
izations excellently agree — the continuum levels agree to bet- 



ter than 10% — except for the OH + line, where the V polariza- 
tion shows ripples and has therefore not been used in the present 
analysis. We combined the data from the three LO tunings and 
the two polarizations to obtain the average spectra shown in Fig. 
[TJ For double-side band mixers with equal sideband gains, fully 
saturated lines are expected to absorb only half of the continuum 
signal. This expectation is borne out in the present data. While 
the OH + line is clearly saturated, both H20 + lines show moderate 
absorption, while the H30 + line only shows weak absorption. 



4. Results 

To understand the structure of the absorption profiles and derive 
column densities of the reactive ions, we normalized the spec- 
tra by the SSB continuum intensities. We constructed empirical 
models of the line profiles by combining multiple velocity com- 
ponents, defined by their central velocity, line width, and opac- 
ity. The predicted spectrum was then computed by combining 
the contributions of the hyperfine components (3 for OH + and 
5 for H20 + as listed in Table [TJ to the overall line profile. We 
show in Fig. [TJ spectra normalized by the continuum intensities, 
with the synthetic spectra superimposed. Table[2]lists the model 
parameters. 

Although we obtained good fits of the OH + and H20 + pro- 
files, it is likely that the solution is not unique because the num- 
ber of velocity components is not well constrained. The OH + 
ground state lines are known with a good accuracy of ~ 1.5 
MHz. Velocities obtained by the empirical model agree well 
with the velocities where deep HCO + (1-0) (Keene et al. 1999, 
Go dard et al 2010) p-H z O and HF (Neufeld et al. 2010) and 
HI dFish e t al. 2003) absorptions are detected. It is remarkable 
that one of the strongest OH + absorption features at Vlsr ~ 7 
km s" 1 coincides with a stron g dip in the HI profi le that has no 
counterpart in CO emission lFalgarone et alJ (120051) . Overall, the 
velocities of the OH + and H20 + velocity components agree well. 
Although the lines are detected in absorption, given the similar- 
ity of their Local Standard of Rest (LSR) velocities with the rest 
velocity of G10. 6-0.4, it is likely that the most blueshifted ve- 
locity components are associated with the G10. 6-0.4 molecular 
complex. 

Below we discuss the column densities and abundances of 
OH + , H20 + , and H30 + in the foreground interstellar gas, with 
LSR velocities higher than 5 kms" 1 . Neufeld et al. (2010) dis- 
cuss the total column density of this material and obtain a col- 
umn density of hydrogen in all forms of Nh ~ 2.7 x 10 22 cm" 2 . 
This figure is consistent with the lower limit derived from the 
oxygen column density measured with ISO, 5 x 10 18 cm" 2 
(Keene et al. 1999), which gives /Vh > 1.7 x 10 22 cm" 2 using 
O/H = 3x 10" 4 . For this line of sight the estimated atomic hydro- 
gen column density is at most /V(HI) = 1.2 X 10 22 cm" 2 (Godard 
et al. 2010), hence the average fraction of hydrogen in molecu- 
lar form is /(i/2) ~ 0.6, with significant variations among the 
velocity components. 

The total OH + and H20 + column densities were obtained by 
summing the contributions of the velocity components with cen- 



troids larger than 7 kms" 1 . We derive /V(OH + ) > 2.5 x 10 14 cm" 2 
and /V(o-H 2 + ) = 4.4 x 10 13 cm" 2 . Although the absorption 
from the p-H20 + line at 607 GHz is blended with emission from 
H 13 CO + (6-5) at 607175 MHz, weak absorption is clearly de- 
tected in the velocity range of the foreground gas along the line 
of sight to G10.6-0.4. The integrated opacity for the LSR veloc- 
ity range 20 to 45 km s" 1 is J Teh - 2kms _1 , corresponding to a 
total column density of ~ 1 x 10 13 cm" 2 . This figure is obtained 
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Fig. 1. HIFI spectra of (a) OH\ (b) o-H 2 + , (c) p-H z O + and 
(d) H30 + . The hyperfine structure of the OH + and o-H20 + lines 
(see Table [TJ is indicated. In both cases, the velocity scale cor- 
responds to the rest frequency of the strongest hyperfine com- 
ponent. In panel c), the velocity scale refers to the frequency of 
the H 13 CO + (6-5) line at 607174.65 MHz. Panel d) uses the fre- 
quency of 9847 1 1 .907 MHz. Panels e) and f) show overlays with 
the synthetic spectra constructed with the empirical models de- 
tailed in Table [2] Grey curves show the profile of the strongest 
hyperfine component. 



Table 2. Empirical models used to analyze the OH + 
H20 + spectra 



and 







OH+ 






H 2 + 




Vlsr 


Av 


T 


N 


A„ 


T 


N 


kms~' 


kms~' 




10 13 cm- 2 


kms -1 




10 12 cm- 2 


-1.0 


14.0 


0.2 


1.3 


14.0 


0.2 


12 


-1.5 


3.0 


0.2 


0.3 


3.0 


0.2 


2.5 


5.5 


6.0 


1.0 


2.8 


3.0 


0.15 


1.9 


8.5 


5.0 


0.8 


1.9 


3 


0.4 


5.0 


11.5 


3.0 


0.5 


0.7 


3.0 


0.3 


3.8 


16.5 


5.0 


1.5 


3.5 








21.5 


5.0 


1.5 


3.5 


10.0 


0.4 


1.7 


29.5 


6.0 


2.0 


5.6 


4.0 


0.8 


1.3 


36.0 


5.0 


2.0 


4.7 


3.5 


0.3 


4.4 


41.0 


6.0 


1.5 


4.2 








45.0 


2.0 


0.5 


0.5 


1.0 


0.1 


0.4 



The OH+ 

hyperfine 



and H 2 + parameters refer to the frequency of the strongest 
component (cf TablefT). 
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Fig. 2. Predictions of the Meudon-PDR model for a slab of A,, = 
1 mag illuminated from both sides. The radiation field and gas 
density are given in each panel. Full lines: n(OH + )/n(H20 + ), dot- 
ted lines: n(OH)/n(OH + ), dashed lines : n(H 3 + )/n(OH + ). Eq. 
(2) determined for each Av in the model is shown in green. 



> 6 and N(H 2 + )/N(H 3 + ) ~ 1.5. These numbers represent the 
average for all velocity components but, given the different ve- 
locity structure of H30 + and OH + , the column density ratio for 
individual velocity components can significantly differ from the 
average value. 



for a limited velocity range, and should therefore be considered a 
lower limit. Given the uncertainties, the data are consistent with 
an o/p ratio of 3 for H20 + , as expected from the spin statis- 
tics. The total column density of H20 + is therefore 7Y(H20 + ) 
~ 6 x 10 13 cm 2 and the abundance ratio 7Y(OH + )/7Y(H 2 + ) > 4. 
Using the total gas column density , we conclude that OH + and 
H 2 + reach fractional abundances, of N(OH + )/Nn > xlO 8 and 
N(H 2 + )/N H ~ 2 x 10~ 9 . 

The H30 + spectrum shows a sharp feature with strongest ab- 
sorption at Vlsr ~ 40 km s and a broad component centered 
at Vlsr ~ 20 kms 1 . This velocity structure is more similar to 
the HF or P-H2O spectra than the OH + and H20 + velocity struc- 
ture, confirming that these molecular ions trace different physi- 
cal conditions. We derived the H30 + column density of the fore- 
ground gas for all velocity components. Assuming an excitation 
temperature of 5 K, appropriate for the conditions of the gas in 
the inner galaxy, we derive 7Y(H30 + ) ~ 4 x 10 13 cirr 2 for the 
LSR velocity range 7-45 kms 1 . Therefore N(OH + )/N(H 3 + ) 



5. Chemistry 

5.1. PDR models 

In order to model the OH + and H20 + chemistry in diffuse 
clouds, we used the Meudon PDR codql a photochemical model 
of a plane parallel and stationary PDR (Le Petit et al. 2006, 
Goicoechea & Le Bourlot 2007). We present models for a rep- 
resentative cloud of Ay=l mag illuminated from both sides. 
Although the total extinction along this line of sight is ~ \Amag 
(Corbel & Eikenberry 2004, Neufeld et al. 2010), the line profile 
clearly shows multiple velocity components, revealing interstel- 
lar clouds along the line of sight. As discussed by Godard et al 
(2010), each of these clouds has low to moderate extinctions, 
corresponding to diffuse or transluscent matter. 

The Meudon PDR code computes the gas temperature gra- 
dient by solving the thermal balance (see Le Petit et al , 2006). 
Our pure gas-phase chemical network is based on a modified 



See http : //pdr . obspm . f r/PDRcode . html 
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Table 3. Main chemical reactions 
Number Reaction 



Rate (cm V) 



OH + + H 2 
H 2 Q + + H 2 



H 2 + + H 
- H 3 Q + + H 



H 2 + + e~ — > products 
H30 + + e~ — > products 
OH + + e~ — » products 



k, = 1.01 x 10 
k 2 = 6.4 x 10 



i-io 



n. 2 — yj.-r /n ±vj 

h = 4.3 x io- 7 (r/300/f)-°- 5 

k 4 = 4.3 x 10- 7 (T/300^)-°' 5 
k s = 3.75 x lQ- 8 (r/300)-° 5 



version of the Ohio State University (OSU) gas-phase network 
updated for photochemical studies. The extinction of UV and 
visible radiation is computed using a dust grain size distribution 
following a power-law (Mathis et al. 1977). The adopted oxygen 
and carbon elemental abundances are 3.02xl(T 4 and 1.38xl0~ 4 , 
respectively. As a starting point we used an ionization rate (of H) 
due to cosmic rays of £# = 2.3xlCT 17 s . 

The oxygen chemistry starts with the 0( 3 Pj) + H + — > + ( 4 S) 
+ H charge transfer process, where reaction rates are sensi- 
tive to the oxygen fine-structure level populations (Stancil et al. 
1999, and references therein). Indeed the charge transfer pro- 
ceeds quickly when the 0( 3 Pi) level (A£'i 2 /A:=228K) is signifi- 
cantly populated, but it can be a slow proceess in cold gas where 
only the 0( 3 P2) ground-state level is populated (Spirko et al. 
2003). As listed in TableEJ OH + is formed by + + H 2 with addi- 
tional contribution from OH photoionization at the cloud edges, 
and by O + H^ . Therefore, the OH + abundance may be sensitive 
to gas temperature, far ultraviolet (FUV) radiation field and £#. 

In order to investigate the possible sources of the high 
OH + /H20 + abundance ratio, Fig. [2] shows photochemical mod- 
els for different densities, and radiation fields (thus different gas 
temperatures). High OH + /H20 + abundance ratios are found in 
the warm gas at low extinction, where the fraction of hydrogen 
in molecular form is significantly smaller than 1 . As a limited set 
reactions need to be considered we built a simplified analytical 
model based on the network listed in Table [3] and find 

«(OH + )/«(H 2 + ) = (k 2 /h) + (k 3 /ki)[n(e-)/n(H 2 )] (1) 

= 0.64 + 430 x (77300r as x [n(e")/n(# 2 )] ■ (2) 

Using f(H 2 ) = 2n(H 2 )/(n(H) + 2n(H 2 )) = 2n(H 2 )/n H ) and as- 
suming that all carbon is ionised (n(e~) - tic - 1.4xlO~ 4 n//)Eq. 
(2)leadston(OH + )/«(H 2 + ) = 0.64+0.12x(77300r°- 5 //(#2)- 
The data suggest that JV(OH + )/JV(H 2 + ) > 4, implying that 
f(H 2 ) < 0.06 at T = 100 K, an appropriate temperature for 
diffuse gas in the molecular ring. Equation (2) determined for 
each Av in the model is shown as green lines in Fig. [2] It fits 
the modeled abundance perfectly, demonstrating that the main 
ingredients have been captured by this simple model. Note that 
the reactions were studied in the laboratory and their rates are 
known with a good accuracy. 

5.2. The role of cosmic rays 

Like H+ (Indriolo et al. 2007), OH + and H 2 + can be used to 
trace the cosmic ray ionization rate for H atoms £#. In the mod- 
erately warm gas where they are produced, the charge exchange 
reaction is rapid. Then, every cosmic ray ionization of H can lead 
to + and we can use the measured OH + abundance to set lim- 
its on £#. In addition to the reaction with H 2 , OH + is destroyed 
in dissociative recombination reactions with electrons. This im- 
plies & > X(OH + )x(e-)n(H)k 5 , ( H > 9.2 x 10~ 12 x n(OH + ) at 
T=100K. 



The lower limit on the OH + abundance in the atomic gas is 
n(OH + )/n(H) > 2 x 10~ 8 , which places a robust lower limit of 
£h > 1.8 x 10~ l9 n(H). This lower limit is very conservative: 
depending on the conditions in the absorbing material, higher 
values of £# may prove necessary to account for the observed 
column densities of OH + and H 2 + . The H^ observations of dif- 
fuse interstellar clouds point toward an enhanced ionization rate 
at low extinction. In a future analysis we will consider how OH + 
and H 2 + can provide additional constraints upon £#. 

5.3. The role of turbulence 

In the warm chemistry driven by intermittent turbulent dissi- 
pation (Godard et al. 2009), the production of OH + and H 2 + 
is enhanced compared to UV-driven chemistry, with an effi- 
ciency depending only on the turbulent dissipation rate, e. For 
e - 2 x 10~ 24 ergcm _3 s _1 , ten times higher than the average 
Galactic value appropriate for the inner Galaxy, the TDR model 
predicts A?(CH + )= 4.3 x 10 13 ctrT 3 , JV(OH + )=2 x 10 12 ctrT 2 
and Af(H 2 + )=2.2 x 10 12 cm~ 2 per magnitude of gas of density 
Yin - 50 ctrT 3 . The ratios between these three species weakly 
depend on e. With the observed CH + column density, 6.7 x 10 14 
ctrT 2 , derived from 13 CH + (Falgarone et al. 2005J3), we there- 
fore expect the following contributions of the TDR chemistry 
in that line of sight: JV(OH + )=3.1 x 10 13 cnr 2 , and A?(H 2 + )= 
3.4 X 10 13 crrT 2 . The OH + column is well below the observed 
value, while that of H 2 + is interestingly close. 

6. Conclusion 

The detection of interstellar OH + , H 2 + , and H30 + provides a 
strong support to our understanding of the gas phase oxygen 
chemistry, in particular the gas phase route to water vapor. The 
combined analysis of the OH + and H 2 + ground state spectra 
has shown that they mostly reside in warmish atomic gas where 
they are produced in a reaction sequence initiated by atomic H 
ionization followed by rapid reaction with traces of H 2 . Their 
abundances and abundance ratio can be used to constrain the 
fraction of hydrogen in molecular form, and to set a lower limit 
for the cosmic ray ionization in the diffuse gas, where few other 
diagnostics are applicable. As these reactive molecular ions can 
be observed from space with the Herschel Space Observatory, 
and in exceptional conditions from the ground, these preliminary 
conclusions should be verified by more extensive observations. 
It is also expected that strong features from both OH + and H 2 + 
will be detected in the spectra of distant galaxies, where elevated 
fluxes of ionizing radiation are produced by active nuclei and in- 
tense star-formation activity. 
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2 In that paper it was underestimated by a factor 8 because of an 
unjustified approximation made in the derivation of the opacity. 
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